Fractalisation drives crystalline states in a frustrated spin system 
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We measure a sequence of quantum Hall-like plateaux at \lq: 9 > q > 2 and plq = 2/9 fractions in the 
magnetisation with increasing magnetic field in the geometrically frustrated spin system SrCu2(BOs)2- We find 
that the entire observed sequence of plateaux is reproduced by solving the Hofstadter problem on the system 
lattice when short-range repulsive interactions are included, thus providing a sterling demonstration of bosons 
confined by a magnetic and lattice potential mimicking fermions in the extreme quantum limit. 

PACS numbers: 75.50.Ee, 75.30.-m, 75.30.Kz, 75.40.Cx, 75.10Jm, 73.43.-f 



Geometrical frustration in the spin dimer material 
SrCu 2 (B0 3 ) 2 JH S S 11 El S H leads to a singlet Shastry- 
Sutherland groundstate at low magnetic fields, but complex 
spin superstructures at higher fields. Our magnetisation mea- 
surements reveal a fine substructure of quantum Hall-like 
plateaux at all 1/q ratios 2 < q < 9 and p/q = 2/9 in 
magnetic fields up to 85 T and temperatures down to 29 mK, 
within the sequence of previously identified plateaux at 1/8, 
1/4, and 1/3 of the saturated magnetization. We identify this 
hierarchy of plateaux as a consequence of confined bosons 
in SrCu2(BC>3)2 mimicking the high magnetic field fractali- 
sation predicted by the Hofstadter butterfly (H for fermionic 
systems. Such an experimental realisation of the Hofstadter 
butterfly has not been previously achieved in real interacting 
materials, given the unachievably high magnetic flux densities 
or large lattice periods required. By a theoretical treatment 
that includes short-range repulsion in the Hofstadter treat- 
ment, stripe-like spin density-modulated phases are revealed 
in SrCu2(BC>3)2 as emergent from a fluidic fractal spectrum. 

The geometrically frustrated spin gap system SrCu2 (603)2 
is unique in its orthogonal arrangement of spin dimers, result- 
ing in an exact direct singlet product incompressible ground- 
state, formally described by Shastry and Sutherland fl. 
Dimers comprising pairs of S = 1/2 spins on neighboring 
Cu 2+ ions bound by an intra-dimer Heisenberg coupling, J, 
are coupled less strongly to orthogonal dimers by an inter- 
dimer coupling J' (Fig. la inset), forming weakly coupled 
layers within the tetragonal crystal structure. Magnetic sus- 
ceptibility and inelastic neutron scattering (INS) experiments 
measure a spin gap of A = 34K, from which J I, J' ~ 1.47- 1.67 
(J = 71.5-100K, J'= 43-68K) [H [H [H is estimated, 
placing the system just within the predicted exact groundstate 
regime J I J 1 > 1.35 13]. In an applied magnetic field H, 
the groundstate incorporates a finite density of spin triplets for 
which the groundstate configuration is no longer known. 

We perform magnetisation measurements on oriented 
SrCu2(B0 3 )2 crystals of ~ 2mm x 2mm x 0.5mm using can- 



tilever magnetometry and ~ 0.3mm x 0.3mm x 2mm using 
pulsed fields, cut from larger single crystals grown using a 
self flux by floating zone image furnace techniques 1U4I1 and 
characterised by INS III 211 and Laue diffraction techniques. 
Measurements up to 35T are made using a 10 thick CuBe 
torque cantilever rotated in-situ to bring the crystalline c-axis 
within 2° of the applied magnetic field, and repeated on an- 
other sample up to 45T using a 50 /im thick cantilever to 
confirm reproducibility and facilitate comparison with pre- 
vious magnetisation measurements ]2|, [16l [171 LL8J]. Torque 
measurements were made in a portable dilution refrigerator 
in continuous magnetic fields at NHMFL, Tallahassee - en- 
abling a level of sensitivity exceeding that possible in pulsed 
magnetic fields. The measured torque is converted to abso- 
lute values of magnetisation by multiplication by a constant 
rescaling factor and a small quadratic background subtraction 
obtained on comparison with measured values of pulsed field 
magnetisation 1161 Il7ll (comparison shown in inset to Fig. 
la). In contrast to ref. Il9ll . we find the magnetisation ob- 
tained from torque measurements to correspond closely to that 
measured by susceptometry in pulsed magnetic fields. Mea- 
surements between 38-85T are made using a coaxially com- 
pensated wire-wound magnetometer during the 10 ms insert 
magnet pulse in the NHMFL (Los Alamos) multishot 100T 
magnet. While the short pulse duration necessarily means 
a relatively large noise floor for these ultra high field mea- 
surements, the identification of prominent features at these 
high fields is facilitated, otherwise unachievable by other tech- 
niques. Pulsed field measurements are made at T=500 mK 
during 2 magnet shots on each of 2 different samples in order 
to confirm reproducibility. 

On applying H within 2° of the crystalline c-axis, a discon- 
tinuous rise in magnetisation m z is observed above a thresh- 
old magnetic field ~ 19.5 T at 29 mK (Fig. la), indicat- 
ing the onset of significant triplet population yj]. On further 
increasing H, a finely-spaced sequence of plateau-like fea- 
tures (these are rounded due to Dzyaloshinskii-Moriya terms 
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FIG. 1: [Colour online] Experimentally measured magnetisation of 
SrCu2(BOs)2 for H || c (a) Extracted from torque measured using a 
lO/im cantilever in static magnetic fields up to 35 T at representative 
temperatures. Plateaux at mjm sa t = plq are indicated by peaks 
in inverse susceptibility AH I AM (red). The lower inset shows that 
low temperature data obtained using a 10/im cantilever up to 35 T 
(green), a 50^im cantilever up to 45 T (pink), ultrashort pulsed field 
measurements from 38 to 85T (blue) and representative published 
pulsed field magnetisation (violet) fl6ll all overlay well. The upper 
inset shows the coupling of Cu 2+ spins within an SrCu2(BO;i)2 layer 
(b) Contour plot of the temperature dependent differential suscepti- 
bility interpolated between measurements using the most sensitive 10 
/an cantilever at 29, 52, 160, 208, 294, 346, 450, 560, 628, 665, 705, 
805 mK, with bright regions indicating transitions between plateaux. 
Dotted lines represent the onset of a significant triplet population, 
and the plateau phase boundary (guide to the eye). Crosses represent 
heat capacity features (in pink) from Ref. |2Ch . and nmr features (in 
cyan) from Ref. 113211 . 



HTot [T5I1 and thermal smearing, we refer to them hereon sim- 
ply as plateaux) appear in m z . Maxima in the inverse differ- 
ential susceptibility correspond to the midpoint of flat plateau 
regions, flanked by thermally rounded transitions. Magnetisa- 
tion plateaux values thus located are found to occur in the se- 
quence m z /m sat = 1/9, 1/8, 1/7, 1/6, 1/5, 2/9 and 1/4 (within 
an error margin of 2%), where m sat is the saturation mag- 
netisation (Fig. la). The corresponding magnetic fields as- 
sociated with each plateau [midpoint (width)] are 27.4(0.4) T, 
28.4(1.2) T, 29.7(0.8) T, 30.6(0.6) T, 32.6(0.8) T, 33.8(0.6) T, 
and 36(6) T respectively - the boundaries being defined by 



maxima in the differential susceptibility. The 1/3 plateau is 
located at 58(24) T from pulsed field measurements (Fig. la 
inset). Fig. lb shows the phase boundary corresponding to 
melting temperatures of the cascade of ordered plateau phases 
located from peaks in the differential susceptibility. The in- 
creased sensitivity of the torque measurements performed in 
continuous magnetic fields and significantly lower tempera- 
tures enables the observance of additional plateaux previously 
unobserved in magnetisation measurements at elevated tem- 
peratures (T > 450 mK for H || c and T > 80 mK for 
H || a) using pulsed magnetic field measurements up to 70T 



171 LUi, which reported only the 1/8, 1/4 and 1/3 



a 

plateaux. Subsequent measurements of nuclear magnetic res- 
onance (nmr) spectra 13211 performed down to temperatures 
of 0.1 9K find differences between the spectral shape at field 
values 27.5 T, 28.7 T and 29.9 T, providing corroboration for 
the fine plateaux hierarchy at 1/9, 1/8 and 1/7 of the satura- 
tion magnetisation identified here by torque measurements. 
While the m z /m sa t= 1/2 plateau was theoretically predicted 
in Ref.Jl l22l 12311 . we observe the first hint of its existence in 
the pulsed field experiment reported here at fields exceeding 
80T. The 1/3 plateau is sufficiently stable and extended in field 
range to enable its unambiguous identification even against a 
relatively large noise floor (due to the short pulse duration in 
the ultra high field 100T magnet), whereas experimental evi- 
dence for the 1/2 plateau is more suggestive in nature. 

We proceed to investigate the origin and character of the 
sequence of field-tuned incompressible spin triplet configu- 
rations associated with the observed magnetisation plateaux. 
The underlying confinement of spin-triplet motion responsible 
for these features arises from the effect of geometrical frustra- 
tion in SrCu2(BC>3)2. The minimal model for describing the 
5 = 1/2 magnetic lattice is the Shastry-Sutherland Hamilto- 
nian: 

U = £ JijBi ■ Sj 9HbHY,S* = £ IfQlbj + b)h) 



(1) 



where i and j denote the lattice sites, Si is the spin 1/2 op- 
erator on site i, and b\, b i are hard-core bosons creation and 
annihilation operators at site i that provide an alternative de- 
scription of a spin system via the Matsubara-Matsuda trans- 
formation 1 24]: Sf = b\ and Sr = 6 ; . The off-site density- 
density interaction described by the last term (n, = b\b i is the 
number operator at site i) corresponds to the Ising term, S* S? , 
in the spin representation. We note that the local magnetiza- 
tion along the field direction (Sf) = (rij) — 1/2 corresponds 
to the particle density in the bosonic representation. 

The computation of m z (H) requires an approach that accu- 
rately captures the underlying hierarchy of energy levels and 
allows to obtain the groundstate energy for a quasi-continuum 
of densities - unlike previously attempted S z = 1 triplet hard- 
core bosonisation 122l 12311 and exact diagonalisation tech- 
niques Ji, 21]. An alternate paradigm is indicated by our ex- 
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perimental observation of plateaux at all 1/q ratios of m sa t 
for 2< q <9 and p/q = 2/9 - reminiscent of the quantum Hall 
effect 1E5U26I1 described by Landau level physics. Hence, we 
begin by adopting a fermionic treatment in which the density- 
density interactions are assumed to be irrelevant ((n.;) = p 
is uniform), as first applied to SrCu2(BC>3)2 by Misguich 
et al. in Ref. 1311. By using a Chern-Simons construc- 
tion on the lattice 112711 . we can map the hard core bosons in 
Eq.O into spinless fermions: b] = f]e l arg(fc ' j)nk , b] = 
g-'Eyj^liilntj. where arg(k,j) is the angle between 
the relative vector, — Yj, and an arbitrary direction. This 
transformation is equivalent to attaching a flux quantum to 
each fermion such that the statistical phase generated from 
fermionic pair-exchange is cancelled by the phase generated 
by the flux quantum via the Aharonov-Bohm effect. After 
this transformation, H becomes a model for gas of spinless 
fermions moving on the same lattice and in the presence of a 
non-local vector potential Ay(r,j). To simplify this problem 
which is rendered intractable due to the non-locality of the 
vector potential, we make the approximation of a uniform flux 
distribution generated by the statistical field H s = 4<J>o^r 01 
(p = J2i( n i) i s th e particle density and N s is the number 
of sites). We thus realise a gas of interacting spinless fermions 
in a strong magnetic field H s . 

The motion of the bosons on the geometrically frustrated 
Shastry-Sutherland lattice is therefore akin to orbital confine- 
ment in two-dimensional quantum Hall systems by a (statisti- 
cal) magnetic field 11251 12611 . Plateaux with a finite Hall con- 
ductance arise when the chemical potential (p = g c pgH) re- 
sides in a gap of the energy spectrum of fermions with aver- 
age density p = m z + 1/2 — a 2 iJ s /4$ D (a 2 is the area of the 
unit cell). The interplay of cyclotron (mean particle separation 
'^2') and lattice 'a' length-scales within this single particle 
band-filling picture was theoretically captured by a Hofstadter 
butterfly of fractalised energy gaps in 1976 Jg]. Minimisation 
of the ground state energy E(m z ) for each value of H yields 
the dependence of m z on H in this non-interacting limit. 
Following this procedure for SrCu2 (603)2 (fractal spectrum 
shown in Fig. 2a), theoretical m z curves are obtained (shown 
in Fig. 2b for J = 70 K, J/J' —2.2 optimised to match ex- 
perimental plateaux lEitl ). the shape of which agrees reason- 
ably well with the measured magnetisation and the values 
of observed plateaux (Fig. 2b) - an advance provided by this 
fermionic treatment over previous models JJS l21, 22, 231. 
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The importance of density-density interactions, however, is 
apparent from the experimental features unexplained by the 
uniform Chern-Simons treatment. While this treatment cap- 
tures most of the measured plateau values and the overall 
shape of m z (H) by neglecting density-density interactions, 
it fails to predict some experimental plateaux. In addition, 
the uniform incompressible liquid groundstate cannot explain 
the broken translational symmetry observed by nmr measure- 
ments at the m z /m sat = 1/8 plateaux J3l or the finite tem- 
perature phase transitions measured using heat capacity I20I1 . 
This suggests that off-site repulsive interactions [see Eq.([T]i] 
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FIG. 2: [Colour online] Comparison between theoretical simulations 
and experiment, (a) Computed fluid Hofstadter spectrum using a 
Chern-Simons treatment of the Shastry-Sutherland lattice in the non- 
interacting limit for J/J'= 2.2 with the Fermi energy in green - 
plateaux in the measured magnetisation correspond to gapped states, 
(a) Simulation of m z lm sa t with J I J' = 2.2 using system sizes ~ 500 
x 500: in the uniform case neglecting interactions (green) and in the 
non-uniform case driven by short-range repulsive interactions (red). 
J / J' was selected to best match the experimental data from a range 
of values incremented by ± 0. 1 . Magnetisation computed from the 
interaction-driven spectrum agrees well with experiment. 



play an important role in modifying groundstate character. To 
include this effect, we consider an unconstrained and self- 
consistent mean field decoupling mrij ~ {n>i)nj + (nj)n>i — 
{rii)(nj), thereby allowing the mean value of the local density, 
(rii), to relax on each site of a super unit cell. Figj3]shows the 
spin density profile corresponding to the lowest energy solu- 
tions of each of the m z lm sat = p/q plateaux, revealing that 
stripe-like density modulations of the spin (except for the 1/2 
plateau) mediated solely by short range interactions signifi- 
cantly lower the energy. 

While the predicted triplet superstructures in Figj3] have 
some features in common with spin configurations predicted 
in Refs.Ji, 22, 23], a qualitative difference is seen in the al- 
ternation of one-dimensional superstructures with arrays of 
spiral-like clusters to accommodate coexisting microscopic 



4 




',V - v > ; v . • »y ' '> r" v'> , - 'V- , V - 



yi\ ■? '.J, 1.7- L£ 



FIG. 3: [Colour online] Minimal energy spin density profiles for 
each plateau. Calculated using several different trial super unit cells 
(sizes ~ (2-8) x (2-36)) with different aspect ratios and on differ- 
ent finite lattices (~ 500x500 original unit cells). Blue to red in- 
tensities represent magnetisation from -1/2 to 1/2. Polarised dimers 
with chiefly singlet contribution (negative magnetisation on one site) 
represented by thin black lines, those with strong triplet contribu- 
tion (positive magnetisation on both sites) represented by thick black 
lines with widths proportional to average magnetisation. A stripe 
super-structure modulates an approximately uniform distribution of 
polarised-dimer clusters for m~ lm aat < 1/3 plateaux. 



length-scales of the lattice 'a' and the mean particle sepa- 



ration (analogous to the cyclotron radius) 



This com- 



mensuration of lengthscales is characteristic to the Hofstadter 
fractal spectrum, the gap hierarchy of which is inextricably 
linked to the distinctive pattern of spin density modulation 
at each incompressible plateau state. Starting with the flu- 
idic spectrum (Fig. 2a), the gap structure is rearranged due 
to a change in lattice potential as each p/q incompressible 
plateau state stabilises triplet stripes separated by an average 
distance of q lattice parameters (where p/q is the density of 
S z = 1 triplets) 12911 . The role of interactions is further am- 
plified by the quasi-degenerate energy level sub-spectrum at 
lower m z lm sat plateaux. For the associated low p/q ratios 
with q >>1 (where q determines the number of sub bands), 
interactions stabilise a unique hierarchy of gaps from the fine 
gap substructure within a bounded spectrum (seen in Fig. 2a). 

Spectacularly, the non-uniform solution corresponding to 
plateaux spin-density profiles in Fig J3] captures the entire se- 
quence of measured magnetisation plateaux, including those 
observed at p/q = 111 and 2/9 which were not predicted by 
the uniform model - comparison with experiment is shown 
in Fig. 2b for a modified interaction-driven spectrum com- 



puted for values of J = 75.5 K and an optimised 'effective' 
ratio J I J' =2.2 that best fit experimental data. The remark- 
able correspondence of predicted and experimental plateaux 
as well as the abrupt step in magnetisation preceding the on- 
set of the plateaux at p/q = 1/9, and the background shape of 
the experimentally measured magnetisation reveals a marked 
improvement over the uniform treatment. 

We thus interpret the observed plateau states in 
SrCu2 (BC>3)2 to represent the emergence of a new type 
of stripe-like crystalline state resulting from the instability of 
the incompressible Hall fluid to interactions. Unlike Wigner 
crystalline states, the structure of these density-modulated 
states is seen to be determined by the density-density sus- 
ceptibility of the uniform solution, as inferred from the p/q 
fractions that are stabilised as a function of field - similar to 
the sequence which appears in the fluid Hofstadter butterfly. 
We are now provided with a new perspective from which to 
question the continued stability of groundstates in metallic 
systems (composed of charged fermions) in the formidably 
high equivalent experimental magnetic field limit. The 
intermediate interaction regime in which the geometrically 
frustrated SrCu 2 (603)2 system lies raises the intriguing 
possibility of a dual groundstate in which crystalline density- 
modulation retains topological properties of a Hall fluid, 
enabling analogies to be made with 'Hall crystals' J30I1 and 
field-induced charge-density wave states I3H1 . 
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